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The beta decay of 192 ' 190 pb has been studied using the total absorption technique at the 
ISOLDE(CERN) facility. The beta-decay strength deduced from the measurements, combined with 
QRPA theoretical calculations, allow us to infer that the ground states of the 192 ' 190 pb isotopes are 
spherical. These results represent the first application of the shape determination method using the 
total absorption technique for heavy nuclei and in a region where there is considerable interest in 
nuclear shapes and shape effects. 


I. INTRODUCTION 

One peculiarity of nuclear systems is that in general 
they are not spherical and take up quite a variety of dif¬ 
ferent shapes. They may have axially deformed prolate 
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and oblate, triaxial, and octupole shapes As a re¬ 

sult the description of nuclear shapes presents a major 
challenge to nuclear models. Here we are concerned with 
shapes in heavy nuclei and in particular with the deter¬ 
mination of shapes using total absorption spectroscopy. 
The method, which is based on a proper measurement 
of the beta strength as a function of excitation energy in 
the daughter nucleus and its comparison with theoretical 
calculations, has only been applied up to now to medium- 
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mass nuclei in the A ~ 80 region SH3. The present result 
illustrates the more general scope of the total absorption 
technique as a tool for inferring the shapes of exotic nuclei 
when combined with the results of shape-sensitive theo¬ 
retical calculations. In this study we have selected the 
neutron-deficient Pb nuclei. They have been the subject 
of intense experimental and theoretical interest in recent 
years. The main reason is the existence of at least one 
0+ excited state below 1 MeV in each of the even-even 
Pb isotopes between A=184 and A=194. This offers ex¬ 
cellent opportunities to study shape effects and shape 
co-existence Q. The cases that have attracted most in¬ 
terest are the 186 ’ 188 pb nuclei where there are even two 
0+ excited states below 700 keV [fll. Theoretically, 
the existence of the low-lying excited 0 + states has been 
interpreted as the result of the combined effect of the 
proton shell gap at Z=82 and the influence of a large 
number of neutron holes below N=126. The role of the 
magic proton number Z=82 is considered to be determi¬ 
nant in the structure of these nuclei. All of them, down 
to 182 Pb are expected to be spherical in their ground 
states HH. 

Theoretically, the co-existence of the low-lying excited 
0+ states has been studied in the framework of a variety 
of models 0. In a shell model picture, the excited 0 + 
states are interpreted as being due to two-quasiparticle 
and four-quasiparticle configurations [13]. On the other 
hand calculations based on phenomenological mean field 
models and the Strutinsky method predict the existence 
of several competin g m inima in the deformation surfaces 
of these nuclei 14, [l5j. Self-consistent mean field cal¬ 
culations [HI EJ and calculations including correlations 
beyond the mean held [3^23| confirm these results. This 
problem has also been studied in the framework of the 
interacting boson model [22], [23} . 

The basis of the present work lies in the fact that it has 
been shown theoretically that the beta decay properties 
of unstable nuclei may depend on the shape of the decay¬ 
ing parent nucleus [24l - [2a ] . In particular cases the calcu¬ 
lated Gamow-Teller strength distribution shows different 
patterns depending on the shape of the parent nucleus 
|26 h 28| . This property can be used to infer the shape of 
the ground state of the decaying nucleus if the Gamow- 
Teller strength distribution has been measured properly. 
For this type of experiment the total absorption spec¬ 
troscopy (TAS) technique is needed, which is the only 
technique that can provide beta intensity distributions 
free from the pandemonium effect 29ij. Measurements 
performed at ISOLDE of the beta decays of the neutron 
deficient 74 Kr and 76 ' 78 Sr nuclei using the total absorp¬ 
tion spectrometer Lucrecia have shown the potential of 
this method mi for nuclei in the A ~ 80 region. 

Theoretical calculations [30, [31| predict that the 
Gamow-Teller strength distributions in the decay of the 
even-even 184 d86,i88,i90,i92,i94pp j so t 0 p es p av e quite dif¬ 
fering patterns depending on the deformation assumed 
for the parent state. This feature can be used to obtain 
new insight into the shape co-existence phenomenon in 


this region if the theoretical results are combined with ac¬ 
curate measurements of the B(GT) distributions in these 
nuclei. This was the primary motivation for the present 
work. In addition, precise measurements of the B(GT) 
distributions in these nuclei are also important to test nu¬ 
clear models further in the Z=82 region. The decays of 
I90,i92pb are 0 f particular interest in this context, since 
the degree of shape mixing in the parent ground state is 
expected to be small |22| . 


II. THE EXPERIMENT 

The measurements were performed at ISOLDE 
(CERN). A 50 gem -2 UCa, target was bombarded with 
a beam of 1.4 GeV protons from the Proton Synchrotron 
Booster (PSB) to produce the radioactive isotopes of in¬ 
terest. The selective laser ionisation source, RILIS [32j], 
was used to enhance the ionisation of Pb isotopes com¬ 
pared to their T1 isobars, which are surface ionised. Then 
the General Purpose Separator (GPS) was employed to 
achieve separation in mass (A). Two measurements were 
made for each mass: one with laser on, to produce a ra¬ 
dioactive beam of the lead isotope of interest, and one 
with laser off, to produce a beam dominated by T1 to 
study possible contaminantion in the Pb measurements. 
The purities of the 192 Pb and 190 Pb beams with the laser- 
on were determined to be 96 % and 93 % respectively. 
In this work we present the results of the analysis of the 
laser-on data for the beta decays of the 192 Pb and 190 Pb 
isotopes. 

Once produced and separated, the beam of the ions 
with the selected mass number was transported to the 
total absorption spectrometer measuring station Lucre¬ 
cia. The beam was implanted in a magnetic tape outside 
the setup, that was then moved in a single step to the 
centre of the TAS detector for the measurements. One 
measuring cycle was composed of symmetrical collection 
and measuring times, which were adjusted to reduce the 
impact of daugther decay contamination in the measure¬ 
ments (cycle lengths: 4 min. for 192 Pb and 2 min. for 
190 Pb), after which the tape was moved again. The pro¬ 
cess was repeated until sufficient data were recorded. 

The TAS detector at ISOLDE, Lucrecia , is made of a 
cylindrically shaped Nal(Tl) mono-crystal with <f> = h = 
38 cm (see the schematic figure of the setup in Ref. 0 )- 
The detector has a transverse hole of <j> = 7.5 cm per¬ 
pendicular to its symmetry axis. The beam pipe is in¬ 
serted in this hole and it is positioned so that the source 
is carried to the centre of the detector for the measure¬ 
ments. On the side opposite to the beampipe, ancillary 
detectors can be placed for the measurement of X rays, 

7 rays, betas, etc. in coincidence with the TAS detec¬ 
tor. The scintillation light from Lucrecia is collected by 

8 photomultipliers. The total efficiency of the detector is 
around 90 % for mono-energetic gamma-rays from 300 to 
3000 keV. During the measurements the counting rate in 
Lucrecia was kept below 10 kHz to reduce the effect of 


3 


second and higher order pileup contributions. Lucrecia 
had an energy resolution of 7.5 % at the 662 keV 137 Cs 
peak. 

In the measurements a germanium telescope was used 
as an ancillary detector. The telescope is composed of 
a Ge planar and a Ge coaxial detector, with the planar 
detector positioned closer to the tape. The planar de¬ 
tector was used with two energy gains: one optimised 
for the measurement of the X-rays emitted in the elec¬ 
tron capture (EC) decay and the other for detecting low- 
energy gamma rays. The coaxial detector was used for 
the measurement of the individual gammas emitted in the 
de-excitation cascades and for looking at gamma-gamma 
coincidences with the planar detector. 


III. RESULTS 

For the analysis of the two decays presented below, a 
TAS gated spectrum was generated by requesting coinci¬ 
dences between the TAS and the T1 X-rays in the planar 
detector (EC component of the Pb decay). Comparison 
of the gates set for different X-ray lines of T1 and Hg 
showed different patterns. This shows the cleanliness of 
the generated X-ray gated TAS spectra (see FigsQ] and 
In addition, the j3 + component of the beta decay 
in the X-ray gated spectra was estimated to be negli¬ 
gible. Accordingly its presence in the X-ray gates was 
discounted. 

The X-ray gated TAS spectra were generated by 
putting gates on the X-ray peaks and subtracting the 
background using gates set on the continuous background 
spectrum as shown in FigsQ] Thus the background and 
pileup, which distort the spectra were taken into account 
in the generation of the spectra for analysis (see Fig[2]for 
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FIG. 1. X-ray gates used for the generation of the analysed 
spectrum for 192 Pb decay. The gates are highlighted in grey 
color. Gates marked with ”Bg” represent the background 
gates used for subtraction. 


the 192 case). 

To analyse the data we need to solve the TAS inverse 
problem: 

jmax 

<k = £ r i 3 ( B )fj (1) 

1=0 

where d, is the content of bin i in the contaminant-free 
TAS spectrum, Rij is the response matrix of the TAS 
setup and represents the probability that a decay that 
feeds level j in the daughter nucleus gives a count in bin 
i of the TAS spectrum, and fj is the beta feeding to the 
level j. The response matrix R^ depends on the TAS 
setup and on the assumed level scheme of the daughter 
nucleus ( B ). To calculate the response matrix first the 
branching ratio matrix B for the levels in the daughter 
nucleus has to be determined. For that the level scheme 
of the daughter nucleus is divided into two regions, a low 
excitation part and a high excitation part. Convention¬ 
ally the levels of the low excitation part and their branch¬ 
ings are taken from high resolution measurements avail¬ 
able in the literature, since it is assumed that the gamma 
branching ratios of these levels are well known. Above a 
certain energy, the cut energy, a continuum of possible 
levels divided in 40 keV bins is assumed. From this en¬ 
ergy up to the Qec value, the statistical model is used to 
generate a branching ratio matrix for the high excitation 
part of the level scheme. The statistical model is based 
on a Back-Shifted Fermi Gas Model level density func¬ 
tion [33| and gamma strength functions of El, Ml, and 
E2 character. The parameters for the gamma strength 
function were taken from [34j and the parameters of the 
level density function were obtained from fits to the data 
available in fui - tdil . As part of the optimisation proce¬ 
dure in the analysis, the cut energy and the parameters of 
the statistical model can be changed. Once the branching 
ratio matrix is defined, the Rij can be calculated recur- 
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FIG. 2. Contributions of the different spectra for the analysed 
spectrum of 192 Pb decay. See Fig. 1 for details on the gates. 
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sively from previously calculated responses using Monte 
Carlo simulations l37H39f . The Monte Carlo simulations 
were validated with measurements of the spectra of well 
known radioactive sources ( 24 Na, 60 Co, 137 Cs). Once the 
R response matrix is obtained, the Expectation Max¬ 
imisation (EM) algorithm is applied to extract the beta 
feeding distributions HO, El| ■ 

The decay level scheme taken from high-resolution 
data was found to be incomplete for both of the cases 
presented in this work. The level scheme of 192 T1 had 
to be modified in order to reproduce the measured TAS 
spectrum. Two levels were added at 494 keV and 690 
keV, based on two gammas seen in the germanium 
detector and on two unexplained peaks appearing 
in the X-ray gated TAS spectrum at those energies. 
When these levels were added to the level scheme, the 
measured spectrum was correctly reproduced in the 
analysis. In this case the cut energy for the statistical 
model was placed at 800 keV in the calculation of the 
branching ratio matrix. A comparison of the analysed 
TAS spectrum with that reconstructed after the analysis 
is presented in Fig. [3] For the analysis of the decay of 
19H Pb the cut energy for the application of the statistical 
model in the level scheme of 190 T1 was placed at 520 keV 
after several trials for this parameter. With this value, a 
similar quality to the 192 case was obtained in the repro¬ 
duction of the analysed spectrum as can be seen in FigU 


The quantity of interest for this study is the beta 
strength. The experimental strength Sp was calculated 
from the beta intensities per energy bin Ip according to 

Eq. H 


Sp(E) 


HE) 

KQeC - E)ti /2 


( 2 ) 


where E is the excitation energy in the daughter nu¬ 
cleus, Qec and ti / 2 are the Q value and half-life of the 
decay and J(Qec — E) is the statistical rate function. 
Our analysis based on the X-ray gated TAS spectra pro¬ 
vides the EC component of the decay. The Ip ( Iec+ Ip+) 
was calculated using the Iec/Ip+ tables of Ref. [42| . 

In the strength calculation for the beta decay of 192 Pb 
into 192 T1 the Qec value of 3310(30) keV |43J and t \/ 2 
= 3.5(1) min [H] were used. The EC+ /3 + branch for 
this decay is 99.9941(7) % [H, [45j (the rest is due to 
alpha decay). For the beta decay of 190 Pb into 190 T1, 
the values of Qec = 3960(50) [t| and ti / 2 = 71(1) s 
Eg} were employed. In this case the EC+ /3 + branch is 
99.60(4) [43, |48| and as in the 192 Pb case, the remainder 
of the decay goes is by alpha emission. In Tables HI and HT1 
the results for the beta feedings from the measurements 
and the deduced beta strength are shown. 


IV. DISCUSSION AND CONCLUSIONS 


The primary aim of this work was to compare the beta 
strength distribution obtained from the TAS measure¬ 
ments with the calculations of Refs. Ham , in order 
to infer the shape of the parent nucleus. The TAS re¬ 
sults indicate that the high resolution studies of these 
decays suffered from the pandemonium effect [29}. The 
measurements showed a redistribution of the beta de¬ 
cay probability compared to earlier high resolution mea¬ 
surements, and a sizable amount of feeding up to the 
Qec value, which is characteristic of this systematic er¬ 
ror. More specifically, the last level observed in high 
resolution studies of the decay of 192 Pb was at 1195 keV 
and above that energy we have detected 41 % feeding. 
Similarly in 192 Pb the last level observed in high resolu¬ 
tion was at 1854 keV and we have detected 18 % feeding 
from that level up to the Qec value. In general terms, 
both the experimental feeding and strength distributions 
of the two decays show a prominent peak at low energy 
(around 1 MeV), which dominated the high resolution 
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FIG. 3. Comparison of the analyzed and reconstructed spec¬ 
trum after the analysis for the 192 Pb decay (a) and the de¬ 
duced feeding distribution (b). 
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TABLE I. Deduced Iec and Iec+p+ values (Ip values) for the levels populated in the decay of 192 Pb for the accepted analysis. 
The last row of the table (total) includes the sum for all the levels up to the last accepted level in the analysis. Sums for two 
energy intervals are also given for the TAS results, one from the cut energy up to the last accepted level in the analysis and the 
other from the last known level from high resolution up to the last accepted level. The high-resolution (HR) data were taken 
from [dil,[45] and the Qec value from [431] . 


Elev 

r 

Iec [%] 


Iec+p+ 

[%] 

B(GT) [g\/ 4tt] 

[keV] 


TAS 

HR 

TAS 

HR 

TAS 

167.49 

i(“) 

~0 

13(9) 

-0 

15(10) 


371.05 

i(-) 


5.8(14) 

-0 

6.3(15) 

-0 

413.98 

(l-,2") 

~0 

- 

-0 

- 

~0 

494^ 


1.4(8) 

- 

1.5(9) 

- 

0.001(1) 

690 ^ 


3.2(22) 

- 

3.3(2.4) 

- 

0.003(2) 

775.67 

(0-.1-) 

8.1(32) 

22(3) 

8.4(33) 

23(3) 

0.007(5) 

1195.46 

1+ 

38(4) 

57(7) 

38(4) 

58(7) 

0.05(2) 

800-3200 


86(5) 

- 

87(6) 

- 

0.24(4) 

1240-3200 


41(2) 

- 

41(2) 

- 

0.18(3) 

TOTAL 


98(7) 

98(9) 

100(7) 

102(1) 

0.24(4) 

a Levels added to the known level scheme to reproduce the experimental TAS spectrum. 



TABLE II. Same 

as Table [I] but for the 190 Pb decay case. The HR data 

2 

were taken from l47l. [48] and the Qec value from 0] • 

El ev 

,r 

Iec [%] 


Iec+p+ 

[%] 

B(G T) [g 2 A /Hr\ 

[keV] 


TAS 

HR 

TAS 

HR 

TAS 

151.31 2. 

(i)- 

~ 0 

< 17 

~ 0 

< 21 

~ 0 

158.15^ 

(0,1)' 


2.1(12) 


2.6(15) 


210.55 

(i)- 

~ 0 

< 17 

- 0 

< 21 

~ 0 

274.17 

(1.2)" 

~ 0 

* b 

- 0 

* b 

~ 0 

372. 75^ 

(O’, I") 


1.3(5" 


1.6(6" 


376.26^ 

(I)" 

~ 0 

4.6(16) 

- 0 

5.4(19) 

~ 0 

416.68 

- 

0.6(3) 

< 19 

0.7(4) 

< 22 

0.0008(4) 

495.07^ 

(I") 


< 19 


< 22 


539.811 

(0,1) 

1.1(5) 

1.7(4) 

1.3(5) 

2.0(4) 

0.002(1) 

598.33 

(1-) 

~ 0 

< 9.7 

- 0 

< 11 

- 0 

738.99 

(0-, 1) 

~ 0 

< 9.9 

- 0 

< 11 

~ 0 

890.72.1 

(0-, 1) 


3.8(6) 


4.1(6) 


942.211. 

1+ 

43.9(10) 

38(5) 

48.5(11) 

41(5) 

0.081(4) 

1235.5 

1 

10.6(3) 

4.3(6) 

11.3(4) 

4.5(6) 

0.024(2) 

1854.5 

1 

5.9(9) 

0.7(2) 

6.0(12) 

0.7(2) 

0.022(2) 

520-3700 

- 

92(1) 

- 

98(1) 

- 

0.32(4) 

1880-3700 

- 

18(1) 

- 

18(1) 

- 

0.14(3) 

TOTAL 


94(2) 

99.6(4) 

100(2) 

100.0(4) 

0.32(4) 


a Levels not resolved by the TAS, their Ip is assigned to levels 151, 376, 540 and 942 respectively. 
b No information available 


data. 

For the theoretical description the formalism described 
in S E| was employed. It uses the quasiparticle- 
random-phase approximation (QRPA), with an axially 
deformed Hartree-Fock (HF) mean field generated by 
Skyrme forces, and pairing treated in the BCS frame¬ 
work. The equilibrium deformation is obtained self- 
consistently as the shape that minimises the energy. The 
potential energy curves obtained with this formalism ex¬ 
hibit oblate, spherical and prolate minima lying close in 
energy that are identified as the ground state and low- 
lying intruder states. Then the HF+BCS+QRPA ap¬ 


proach is applied to the calculation of the B(GT) dis¬ 
tribution, assuming that the parent and daughter nuclei 
have the same deformation. The B(GT) distributions 
of the 184 ~ 1!,4 pb isotopes clearly show different patterns 
depending on the assumed deformation of the parent 
ground state ( Fig. 5 of ([305] ), a feature that can be used 
for shape determination if combined with a proper de¬ 
termination of the beta strength distribution. It is worth 
noting that while the relative equilibrium energies are 
sensitive to both Skyrme and pairing interactions [30j |. 
the B(GT) distributions are more sensitive to the nu¬ 
clear deformation than to those interactions j3(|. In the 
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following discussion the theoretical results of [UJ using a 
Sly4 Skyrme force are used. 

Although lead nuclei are expected to be spherical in 
their ground states [TTj | . the use of this method can be 
viewed as a test to confirm its validity in this region, 
with the idea of extending it to other possible cases of 
interest such as Pt, Po and Hg, which are not necessarily 
spherical in their beta-decaying states. We remind the 
reader, that this method has only been applied so far 
in the A ~ 80 region. With the availability of a reliable 
technique to measure the B(GT) strength over the whole 
Qec window, the possibility is open also for further tests 
of nuclear models in the region around Z = 82. 

One of the best ways to compare theory and experi¬ 
ment is using the accumulated beta strength. The accu¬ 
mulated strength at each energy was obtained by adding 
the strength observed up to that excitation energy. The 
results are presented in Fig. [5] In the comparison the 
theoretical results from Ref. 3l| have been scaled by a 
standard quenching factor (9a/ 9v)eff = ^>-77 (g a/ 9v)i 
that is, a global factor of 0.77 2 ~ 0.6. The theoret¬ 
ical calculations produce sizeable fragmentation of the 
strength only in the deformed cases, which is reflected 
by the “staircase” character of the accumulated strength. 
The two panels show a nice agreement with the theoret¬ 
ical results (total strength) for spherical ground states, 
but showing less spreading or fragmentation than in the 
A ~ 80 region cases a feature partially related to 

the spherical character of the nuclei studied. In partic¬ 
ular the total B(GT) strength to levels below the Qec 
window is very well reproduced. 

In the spherical case, if we compare the theoretical 
and experimental excitation energies where the strongest 
populated peak appears, we see an energy difference of 
approximately 0.7 MeV that is mostly related to the dif¬ 
ficulties of these forces in reproducing single-particle en¬ 
ergies accurately. Skyrme interactions were designed to 
provide universal descriptions of nuclear properties over 
the whole nuclear chart and detailed spectroscopy is be¬ 
yond the scope of the present approach. When moving 
to more exotic neutron-deficient Pb isotopes the mea¬ 
surements indicate that the energy of the state where 
most of the strength is concentrated decreases and the 
total strength increases (including the 188 Pb decay case 
not presented here [II]). These general trends are well 
accounted for by the calculations. 

The quantitative description of mixed shapes is be¬ 
yond the mean field approach used at present, and one 
should consider that the experimental results are in nice 
agreement with theory because of the expected small de¬ 
gree of mixing. Based on the present results we can in¬ 
fer spherical character for the 190 > 192 p& ground states. 
This represents an additional and independent proof of 
the results obtained in [llj| . where small deviations from 


the liquid drop model for these nuclei were interpreted 
as the effect of tiny admixtures of intruder contributions 
in the ground state wave functions, but essentially they 
were considered to be spherical. Another conclusion of 
our work is that the method can be applied to heavy nu¬ 
clei and suggests that this technique (TAS measurements 
combined with theoretical calculations) is suitable for use 
with other deformed nuclei in the region such as the Hg 
and Pt isotopes, which are of particular interest because 
of the existence of isomeric states with different shapes. 

This work has been supported by the Span¬ 
ish MINECO through projects FPA2008-06419-C02-01, 
FPA2011-24553, FPA2013-41267-P, FPA2014-52823-C2- 
1-P, FRPA2012-32443, FIS2011-23565 and CSD-2007- 
00042 (CPAN Consolider, Ingenio2010), by the Euro¬ 
pean Union Seventh Framework via ENSAR (contract 
no. 262010), contract PIOF-GA-2011-298364 and by 
OTKA contracts K100835 and K106035. 
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FIG. 4. (a) Comparison of the analyzed and reconstructed 
spectrum after the analysis for the 190 Pb decay and (b) the 
feeding distribution deduced from the analysis. 
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FIG. 5. Comparison of the accumulated strength for the beta 
decays of 192 ’ 190 pft with theory. The theoretical results are 
presented in two ways: the accumulated direct results from 
the calculations and accumulated strengths calculated with 
folded distributions using 1 MeV width Gaussians on the dis¬ 
crete spectrum. 
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